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We derive the quantum trajectory or stochastic (conditional) master equation for a single su- 
perconducting Cooper-pair box (SCB) charge qubit measured by a single-electron transistor (SET) 
detector. This stochastic master equation describes the random evolution of the measured SCB 
qubit density matrix which both conditions and is conditioned on a particular realization of the 
measured electron tunneling events through the SET junctions. Hence it can be regarded as a 
Monte Carlo method that allows us to simulate the continuous quantum measurement process. We 
show that the master equation for the "partially" reduced density matrix [Y. Makhlin et.al., Phys. 
Rev. Lett. 85, 4578 (2000)] can be obtained when a "partial" average is taken on the stochastic 
master equation over the fine grained measurement records of the tunneling events in the SET. Fi- 
nally, we present some Monte Carlo simulation results for the SCB/SET measurement process. We 
also analyze the probability distribution P{m,t) of finding m electrons that have tunneled into the 
drain of the SET in time t to demonstrate the connection between the quantum trajectory approach 
and the "partially" reduced density matrix approach. 
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I. INTRODUCTION 

The single-electron transistor (SET) is a highly charge- 
sensitive electro-meter and has been suggested as a 
readout device for solid-state charge qubitgJ-'- or spin 



qubits^ (through a measurement of a spin-dependent 
charge transfer). The problem of a charge qubit subject 
to a measurement by a SET has been extensively stud- 
ied in Refs. Q and|3- We refer to the approach of these 
papers as the master equation method of the "partially" 
reduced density matrix. In this approach, one takes a 
trace over environmental (detector) microscopic degrees 
of the freedom but keeps track of the number of electrons, 
m{t), that have tunneled through the SET into the drain 
during time t in the "partially" reduced density matrix. 
If experimentally the number of accumulated electrons 
or current passing through the SET is measured, this ap- 
proach can provide us with information about the initial 
qubit state. But, the system dynamics in this approach 
is still deterministic; i.e., this approach is still in an en- 
semble and time average sense. 

A Monte Carlo method^ which allows one to follow 
each electron tunneling event has been successfully ap- 
plied to simulate transport properties of a SET or more 
complicated single electronics circuits. This method gives 
physical insight into the processes taking place in the 
simulated system. But to our knowledge, it has not yet 
been formally applied to quantum measurement prob- 
lems by a SET detector. In this paper, we provide such 
an investigation. We derive the quantum-jump stochas- 
tic master equation (or quantum trajectory equation) for 
a single superconducting Cooper-pair box (SCB) charge 
qubit (generalization to other charge qubit case is sim- 
ple) continuously measured by a SET. This stochastic 
master equation describes the random evolution of the 
measured SCB qubit density matrix which both condi- 



tions and is conditioned on a particular realization of 
the measured electron tunneling events through the SET 
junctions. We can regard it as a Monte Carlo method 
that allows us to simulate the continuous quantum mea- 
surement process of a charge qubit by a SET. This quan- 
tum trajectory approach (or Bayesian formalism) was 
introduced recentlj5§iZii to describe a charge qubit mea- 
sured by a low-transparency point contact detector. Here 
we present the quantum-jump stochastic master equa- 
tion for the SET detector. Especially, we show that 
the master equation for the "partially" reduced density 
matrix (a "partial" course-grain description) presented 
in Refs. Q and y can be obtained by taking a "partial" 
average on the stochastic master equation over the fine 
grained measurement records of the tunneling events in 
the SET. Finally, we present some Monte Carlo simula- 
tion results for the SCB/SET measurement process. We 
also analyze an important ensemble quantity for an initial 
qubit state readout experiment, P{m, t) the probability 
distribution^^ of finding m electrons that have tunneled 
into the drain of the SET in time t. This analysis demon- 
strates further the connection between the quantum tra- 
jectory approach presented here and the "partially" re- 
duced density matrix approach in Refs.|l| andy. 



II. MODEL HAMILTONIAN 

The Hamiltonian of the SCB/SET system is described 
in Refs. Q and|2 as: 

'H = HsET + T^L + Ti-R + "Hi + Ht + 'Hqb + Hint- (1) 
Briefly, 



WsET — £^SEt(^ — Ng 



(2) 



describes the charging energy of the SET. The charge on 
the middle island is eTV, and the induced charge eNg is 
determined by the gate voltage Vg and other voltages in 
the circuit. 
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(3) 



where r = L, R, I, describes microscopic degrees of free- 
dom of noninteracting electrons in the two leads (left and 
right) and the middle island of the SET, respectively. 
To make the charge transfer explicit, two "macroscopic" 
operators, e^"^ and e^*'^ are included in the tunneling 
HamiltonianiiS in the SET: 



"^T = ^ T^k'sctl^i's 



-i<p 



kk' s 



+ ^T,V.cf,U'.e-^"e^'^+H.c. . (4) 



k' k" s 



The effective Hamiltonian of the uncoupled qubit, writ- 
ten in the charge eigen basis of the number n of extra 
Cooper pair on the island of SCB, is 



^qb ^ ■^(-^chO'z ^ Eja^), 



(5) 



where h — (I — ^z)/'^ with eigenvalues n = or 1. The 
capacitive Coulomb coupling between the charge on the 
SET island and that on the SCB qubit is represented by 



Hint = 2Ei,,tNn 



(6) 



We will consider the case that the leading tunneling 
process in the SET are sequential transitions between 
two adjacent charge states N and iV -I- 1 (say, iV = 
and A'' -I- 1 = 1 states to represent the extra charge on 
the SET island). This would be the case if the applied 
transport voltage across the SET is not too high and 
the temperature is low (for simplicity, we consider the 
zero temperature case). Since, effectively, only two ad- 
jacent charge states N — and A^ + 1 = 1 are con- 
sidered, the charge transfer operators e'^^'^, in this case, 
satisfy e-'-^lA^) = = e''^\N + 1), e"^\N) = \N + 1) and 
g-«0|7V -I- 1) = |A^). The other set of charge transfer op- 
erators satisfy e '^\m) = |?7i ± 1), where m represents 
the number of electrons that has tunneled into the right 
lead (drain) of the SET. 



III. MEASUREMENT RECORDS AND 
CONDITIONAL DENSITY MATRIX 

To be able to describe the measured qubit in a pure 
state continuously, one needs to have the maximum 
knowledge about the change of its state. When the qubit 
interacts with (is measured by) the SET, this informa- 
tion is lost to the SET. For example, each time when 
an electron tunnels onto or off the SET island, it will 
cause a change (e.g., a phase shift) of the qubit state. 



One can recover this information lost, provided that a 
detailed measurement record from the SET is available. 
The transport of electrons through the SET occurs via 
real states of the central island, from N ^ N + 1 ^ N. 
The information of detecting the to'^ electron just tun- 
neling into the drain only tells us that the island state 
now is in |A^) state. Thus knowing the "partially" re- 
duced density matrix ppf{m, t) at every timsi;^ does not 
provide us with the full information. 

One can imagine that in the transport process, elec- 
trons may spend different times in the intermediate 
\N + 1) state, causing different phase shifts to the qubit. 
If the record of the times when electrons tunneling onto 
and off the island is not available from the measurement 
results of the SET, our knowledge about the precise qubit 
state decreases. When this happens, averaging the ran- 
dom dwelling times of electrons on the island over a pe- 
riod of time or over an ensemble of systems will then 
lead to the decoherence of the qubit. Hence, one needs 
to have a measurement record which records whether or 
not an electron tunnels onto or off the central island of 
the SET at each time interval dt. This time interval dt 
should be much smaller than the typical qubit system 
evolution or response time so that no information is lost 
as far as the qubit system evolution is concerned. In this 
sense, effectively the qubit is continuously monitored or 
measured. 

For this purpose, we introduce dNLc{t) and dNj^dt) 
to represent, in the quantum-jump case^, the number (ei- 
ther zero or one) of tunneling events seen in infinitesimal 
time dt through the left and right junctions of the SET, 
respectively. Throughout the paper, the subscript or su- 
perscript c indicates that the quantity to which it is at- 
tached is conditioned on previous measurement results. 
If no tunneling electron is detected, the result is null, 
i.e., dNLc{t) — and dNfidt) — 0. If there is detection 
of a tunneling electron in time dt, then dN^dt) — 1 or 
dNfic{t) = 1. We can think of dNjidt) as the increment 
in the number of electrons Nfidt) — J2 '^^Rc{t) passing 
through the right junction of the SET into the drain in 
the infinitesimal time dt. It is the variable Nfi{t) = rn(t), 
the accumulated electron number transmitted through 
the SET in the drain, which is used in Refs. Q and Q- 
Since the nature of detection results is classical and that 
of electrons tunneling through the SET is stochastic, 
dNi^c/Rcit) should represent a classical random process. 
The measurement record in each single run of experiment 

is the set of times {tl } and {i^ } when electrons tun- 
nel onto or off the SET island, respectively [i.e., ones of 
dNLc{t) and dNjidt) over the entire detection time; see, 
e.g., Fig. mi) and (j)]. 

At first, one may expect that at the end of each 
time interval dt, there are four possible measurement 
outcomes, dNLd*) [1 - dNadt)], dNndt) [1 - dNLdt)], 
[1-dNLdt)] [l-dNRdt)], and dNLdt) dNndt)- It is im- 
portant to realize that a null result (e.g., dN^c = dNjic = 
0) in a time interval dt is still a measurement result or 
outcome. Let us consider the case in the sequential tun- 



neling dominated regime that the probabiUty of electrons 
tunnehng onto and off the SET island within the same 
infinitesimal time interval dt is rather small. In fact, 
the respective probability of dNLc{t) or dNfidt) equal 
to unity is proportional to dt [see Eqs. (I27II and H28() ]. 
Thus the product of dN Lc{t) dN ]ic{t) = 1 occurs with 
probability proportional to dt^ . Since we shall keep only 
terms to order dt in the master equations, we can ne- 
glect the case that dNLc{t) and dNjidt) both equal one 
within the same infinitesimal time interval. The possible 
measurement outcomes then become: dN^dt), dNjidt) 
and [1 — dNLc{t) — dNjidt)]. The first two terms, in this 
case, represent that an electron tunneling event through, 
respectively, the left and right junctions of SET is mea- 
sured at the end of the time interval [t,t + dt). While 
the last term [1 — dN^cit) — dNiic{t)] represents that no 
tunneling event is observed in [t, t + dt). Thus the condi- 
tioned density matrix Wc{t + dt) to order dt at the end 
of the time interval [i, t + dt) can be written as 



Wc{t + dt) = dNLcit) 

+dNRdt) 



WLlcjt + dt) 

Tr[WLic{t + dt)] 

WRlc{t + dt) 



TT[WRiAt + dt)] 



+ [1 - dNLcit) - dNRdt)] 



Wacjt + dt) 
Tr[Woc{t + dt)] 



,(7) 



where WliS + dt), WRic{t + dt), and Woc{t + dt) are 
the unnormalized density matrices, given that an elec- 
tron tunneling event through left or right junction of the 
SET island, or no tunneling event is measured at the end 
of the time interval [t, t + dt). Equation |(Zj) simply states 
that when dN^c = 1 and dNRc = 0, the normalized con- 
ditioned density matrix is WLic{t + dt) /TT[WLicit + dt)], 
and so on. Self-consistently, the ensemble averages 
E[dNLc{t)] and E[dNRc{t)] of the classical stochastic 
processes dNLcit) and dNRc{t) should equal respectively 
the probabilities (quantum average) of electrons tunnel- 
ing through the left and right junctions of the SET is- 
land in time dt, i.e., E[dNLc{t)] = Tr[WLic(t + dt)] and 
E[dNRc{t)]^Tr[WRXc{t + dt)]. 

Formally, we can write the currents through the junc- 
tions as 



lLc{t)^e[dNLc{t)/dt], 

lRc{t)^e[dNRcit)/dt]. 



(8) 
(9) 



The question now is to find expressions for WLic{t + dt), 
WRic(t+dt), and Woc{i+dt) in the model. To do this, we 
derive the unconditional master equation and then use it 
to find WLic{t + dt), WRic{t + dt), and Wocit + dt). 



IV. STOCHASTIC MASTER EQUATION 



Following the same assumptions and approximations 
in Refs. Jj and y and similar derivations in Refs. 7 and 



|g|, we first derive the master equation of "partially" re- 
duced density matrix. By tracing out the microscopic 
degrees of freedom of the left and right leads and the 
island of the SET, but keeping the electron transfer op- 
erators explicitly^'^ (so that we can keep track of effects 
of electron tunneling events on the system density ma- 
trix), we obtain the Born-Markov master equation for 
the "partially" reduced density matrix operator 'W{t) of 
the SCB/SET system (consisting of the qubit, and the 
island and drain of the SET) as: 

[dW{t)ldt] = -(l/?i)[Hqb+Hi„t,W(i)] 

+TLV[e">'{l - h)]W{t) + T'LV[e"t'h]W{t) 

+TRV[e-''f'e''^{l ~ n)]-W{t) + r'^2?[e-^'^e^^n]>V(i) 

-{TL+T'L)[n,[n,e^*W{t)e-^f']]/2 

-(rfl, + r'fl) [n, [n, e-'^e^ wit) e-'^e%l2 , (10) 

where T> is defined for arbitrary operators B and W as 

V{B]W = J{B\W-A{B]W, (11) 

J{B\W = BWB\ (12) 

A[B]W = iB'<BW + WB^B)/2. (13) 

The rates Tl/r and ^'l/r represent the tunneling rates 
(in the left or right junction) with and without the pres- 
ence of the extra Cooper pair on the island of the SCB 
(i.e., n = 1 or n = 0), respectively. They are determined 
by the chemical potentials fJ-L/R. of the leads and the in- 
duced charge Ng on the SET's island: 



Tl = 

Tr = 

r'L = 

r'p = 



i2^TaL/n)[^iL-{l-2Ng)EsET], (14) 

i2TraR/n) [(1 - 2Ng)EsET - m], (15) 

Tl - {47:aLE,^t/h), (16) 

Tn + {A7raREint/n), (17) 



where ol/r = RQ/{4:Tr^RL/R), Rq = h/e^ is the resis- 
tance quantum, and Rl/r represents the resistance of 
the left or right junction. 

The master equation (|10|l has a translational symme- 
try in m space. So by summing all possible values of m 
of the right reservoir (drain) states completely, a closed 
form of the master equation of W{t) = X]m("^ 1^(^)1"^) 
can be obtained. The resultant equation is equivalent 
to Eq. H1U|I but with the replacements of e***^ -^ 1 and 
yV(i) -^ W{t). One may expect to apply the similar 
sum procedure to the island states. However, since ef- 
fectively only two extra adjecent charge states |7V) and 
|7V-|- 1) are considered, a closed form of the master equa- 
tion for the qubit density matrix operator alone, p(t) — 
{N\W{t)\N) + {N + l\Wit)\N+l) = pNit) + pN+iit), can 
not be obtained without further approximations, where 
PN/N+iit) is each a 2 x 2 operator in the qubit ba- 
sis. One approach^ is to assume extremely asymmet- 
ric tunnel junctions for the SET, i.e., one of the tunnel- 
ing rates through the junctions is much larger than the 
other. In this case one can apply the adiabatic elim- 
ination procedure^ to eliminate the degrees of freedom 



of the SET island to obtain the reduced density matrix 
for the qubit alone. But this asymmetric assumption 
is equivalent to treating the SET as effectively a sin- 
gle junction device, similar to a point contact detector. 
Here, however, we take the joint density matrix of the 
qubit and extra charge on the SET island as the system 
density matrix in Eq. {TJ. After evaluating Wc{t) [or 
p%i(t) and p'lfjAf)] from the conditional master equa- 
tion [see Eqs. (|21(l and (|^ ]. we can then find the con- 
ditional qubit density matrix operator alone by writing 
p'^it) = TT^[W,{t)] = pUt) + P%+,{t). 

Using the definition of the superoperator V and the 
fact that the charge transfer operators are explicitly kept 
in each term in Eq. (|l()|l . one can then findi from there 
[or more precisely from the master equation for M^(i)] 
the unnormalized density matrices, given that an electron 
tunneling event through left or right junction of the SET 
island takes place at the end of the time interval [t, t+dt), 
as: 

WLic{t + dt) = dt{rLJ[e'^{l-n)]Wcit) 

-(Ei +r'J[n, [n,e''f'Wcit)e-'*]]/2} ,(18) 
Wnicit + dt) - dt{rRj[e-''^{l-n)]Wcit) 

-{rn + T'j,)[n, [n,e-''^W,{t)e''^]]/2} (19) 

It is required that the unconditional (ensemble averaged) 
density matrix ^[M/c(i + di)] = W{t + dt) = WLiit + dt) + 
Wri {t + dt) + Wo{t + dt). Hence we find Woc{t + dt), from 
the master equation for W{t), as 

Woc{t + dt) = Wc{t)-dt{l/h)[Hq^+ninuWa{t)] 

-dt {TLA[e"''{l - n)] + T'^A[e"t'h] 
+ri,^[e-^'^(l - h)] + T'j,A[e-''>'n]]W,{t){2{)) 

Substituting Eqs. iQHll-lEni) into Eq. Q and replacing 
T:i[Woc{t + dt)] ^\-TY[WLic{t + dt)]-Tr[WLic{t + dt)l 
then keeping only the terms to order dt in the resultant 
equationl, and finally evaluating this equation in \N) and 
A^-|-l) states respectively, we obtain the conditional mas- 
ter equation: 

dp%{t + dt) = -[dNLc{t)+dNR,{t)]p%{t) 

+dNR,{t) [tRP%+^{t)/rRUt)] 

-dt{{l/fl)[H^^,p%{t)]+tLp'},{t) 

-[pLic(i)+^mc(i)]p^(i)} , (21) 
dp'k+iit + dt) = -[dNLcit) + dNR,{t)]p%+i{t) 

+dNLc{t) [fLPN{t)IVLlc{t)] 

~dt{{i/h) [T^qb + 2E„un, Pn+1 (01 
+fflp5^+i(t) 

~[rLic{t) + PRUt)VN+At)} ,(22) 

where VLc{t) and VRc{t) appearing in Eqs. H21|l and H22|) 
are due to the normalization requirement for the density 



matrix after each detection interval dt as in Eq. jT)), and 
are given by 

VLic{t) = rLTr[p5v(i)] + (rL - T^)Tv[np%{t)] , (23) 
VRUt) = rRTr[p^+i(i)] + {V^ - rR)Tr[np^+i(^4) 

The rates T l and T r are defined as 

TlpL = rip^ + (ri-r'J{n,p5v}/2, (25) 

TrP^n+i = rRP^+i-l-(ri,-ry{n,p^+J/2.(26) 

Self-consistently, E[dNLc{t)] and E[dNRc{t)] should 
equal their respective quantum averages, and from Eqs. 
(^ and (QSl) can be written a£^° 

E[dNUt)] = T^[WLicit + dt)] = VLic{t)dt , (27) 
E[dNR,,{t)] = TY[WRic{t + dt)] = VRic{t)dt , (28) 

where VLc{t) and VRc{t) are defined in Eqs. (|23|l and 
(|^. Equations l(?T)) - (P5|l are the main results of the 
paper. One can use them to simulate the conditional 
(stochastic) qubit dynamics under continuous quantum 
measurements by the SET. We will present some simula- 
tion results in Sec. EU 



V. CONNECTION TO "PARTIALLY" 
REDUCED DENSITY MATRIX 

Next, we show that the master equation of the "par- 
tially" reduced density matrix, e.g., Eq. (2) of Ref.^, can 
be obtained^ by taking a "partial" average on Eqs. Il21() 
and (|22|l . First, performing a full ensemble average over 
the observed stochastic process on Eqs. (|21|l and H22|l by 
replacing E[dNLc{t)] and E[dNRc{t)] by their expected 
values Eqs. ^7^ and ^^, and setting E[p'j^{t)] = pN{t), 
we obtain the master equation for the reduced density 
matrix pN{t) and pN+i{t) as 



d 
dt 



PN{t) \ ,i ( [T-f-qb, PN{t)] 

PN+l{t) J n\ [Hqb + 2^intn, pN+l{t)] 

-Tl ^R \f PN{t) 
fi -Tr \pN+l{t) 



(29) 



Then to keep track of the number of electrons m — Nr 
that have tunneled into the drain, we need to identify 
the terms in Eq. (|29|) . which come from Eqs. H21|) and 
(|22|l and have effects corresponding to an electron tun- 
neling through the right junction of the SET. Only one 
such term, originating from dNRc{t) [TRP%+iit)/PRc{t)] 
in Eq. (|21|l . survives in Eq. H29|) . It is in the upper right 
corner of the matrix on the right hand side of Eq. H29|l. 
If m electrons have tunneled through the right junction 
of the SET at time t + dt, then the accumulated number 
of electrons in the drain at the earlier time t, due to the 
contribution of the jump term through the right junction, 
should be {m — 1). Hence, after writing out the number 



(30) 



dependence m or (m — 1) explicitly for the density ma- 
trix in Eq. (|29|l . we obtainii the master equation for the 
"partially" reduced density matrix as: 

d_f pN{m,t) \ ^ i_ ( [Hqh,PN{m,t)] 

dt \ PN+i {m,t) J h\ [Hqb + 2-Bintn, PN+i{m, t)] 

^ f -TLPNim, t) + tiiPN+i{m - l,t) 
\ tLPN{m,t) -rHPN+i{m,t) 

Making a Fourier transform PN/N+i{k,t) = 
Sme~*'""P7V/7V+i(TO,*) on Eq. (1201), we find that 
the resultant equation is exactly the same as Eq. (2) of 
Ref. |3- If the sum over all possible values of m is taken on 
the "partially" reduced density matrix [i.e., tracing out 
the detector states completely, pN{t) = X^m Z'^l™'*)]' 
Eq. H3U|I then reduces to the master equation of the 
reduced density matrix, Eq. (|29|l . This procedure of 
reducing Eqs. ^^ and ^^ to Eq. ^^ and then to 
Eq. (|29|l . by successively disregarding information that 
distinguishes different states of the detector, provides a 
connection between the approach of Refs. [l| and and 
the more detailed stochastic master equation used here. 
To further demonstrate this connection, we analyze in 
the next section an important ensemble quantity for 
an initial qubit state readout experiment, P{m,t), the 
probability distribution of finding m(= Nfi) electrons 
that have tunneled into the drain in time t, considered 
in Ref. il| 



VI. SIMULATION RESULTS AND 
DISCUSSIONS 

Although the "partially" reduced density matrix 
approacbii^ can provide the information about the ini- 
tial qubit state, the system dynamics in this approach is 
still deterministic; i.e., this approach is still in an ensem- 
ble and time average sense. If a measurement is made 
on the qubit system and the results are available, the 
state or density matrix is conditioned on the measure- 
ment results. If the subsequent system evolution after the 
measurement is concerned, the conditional or quantum 
trajectory approach should be employed. In particular, 
to describe the conditional dynamics of the qubit sys- 
tem in a single realization of continuous measurements, 
which reflects the stochastic nature of electrons tunnel- 
ing through the SET junctions, we should use the condi- 
tional, stochastic master equations H21|) - (|28|l . 

A set of typical quantum trajectories, generated us- 
ing Eqs. lf^ - (PS|) . is shown in Fig.^a)-(h) and its cor- 
responding randomly distributed moments of detections 
are presented in (i) and (j), where p"" = {n\p\n'). Due 
to Coulomb blockade, when an electron is on the SET 
island, Eqs. (|21|I - H28I) ensure that no electron can tunnel 
onto the SET island, i.e., guarantee dN^c = for the next 
time interval dt. Note that in this case, we still have two 
possible measurement outcomes of either dN^c = 1 or 
in the next time interval dt. But only after a detection of 




FIG. 1: Set of typical quantum trajectories and correspond- 
ing detection record for the initial conditions: qubit state in 
(V^jO) -I- |l))/2 and the SET island in iV = state. Other 
parameters are J5int = l.SBch = ISOOBj, Fl = (3/4)Fi{ = 
3-Bch/ft, OIL ~ ciR = 0.03 and the time is in unit of Ti/ Ech- 



dNjic = 1 in the next or the next several time interval(s), 
could another electron tunnel onto the SET island again. 
This is the measurement record shown in Fig.^i) and (j). 
They are in the order of exactly alternating dN^c = 1 
and dNjic — 1 time sequence. The conditional evolutions 
of the qubit alone shown in Fig. ^e) and (f) can be ob- 
tained from the sum of the joint state evolutions of (a) 
and (b), or (c) and (d) , respectively. The probabilities. 



Pr 



0/1, c 



Trqb[po/i,c 
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pj}^^^, of the SET island 



state alone in (g) and (h) can be obtained by summing 
the evolutions in (a) and (c), or (b) and (d), respectively. 
The conditional evolutions in (a)-(h) differ considerably 
from their ensemble average counterparts. 

In this conditional or quantum trajectory approach, we 
are propagating in parallel the information of the condi- 
tioned (stochastic) state evolution and detection record 
of dNLc{t) and dNRc{t) in a single run of a continuous 
measurement process. One can see that in this approach 
the instantaneous system state conditions the measured 
electron tunneling events through the SET junctions [see 
Eqs. (|27|) and (|28|) ]. while the measured electron tun- 
neling events themselves condition the future evolution 
of the measured system [see Eqs. (|21|) and (|22() ] in a 
self-consistent manner. Each set of quantum trajectories 
(stochastic state evolutions), obtained from the stochas- 
tic master equations (|21f) - (l28|l . mimics a single history of 
the qubit state in a single run of the continuous measure- 
ment experiment. The stochastic element in the quan- 
tum trajectories corresponds exactly to the consequence 
of the random outcomes of the detection record of the 
tunneling events in the SET. The macroscopic ensem- 
ble measurement properties can be calculated by using 
large ensembles of single electron tunneling events (fine 
grained measurement records). 
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FIG. 2; Probability distribution P{m,t = 2000). The initial 
conditions and parameters are the same as those in Fig. Q 



If only one measurement value is recorded in each run 
of experiments [for example, the number of electrons m 
that have tunneled into the right lead (drain) in time t] 
and ensemble average properties [for example, P(m,t), 
the probability distribution of finding m = (Nn) elec- 
trons that have tunneled through the right junction into 
the drain in time t] are studiedi*^ over many repeated ex- 
periments, the quantum trajectory approach will give the 
same result as the master equation approach of the "par- 
tially" reduced density matris^. However, more physical 
insight can be gained in the approach of quantum trajec- 
tories. We demonstrate this feature below. 

We consider the case that the SCB/SET system is in 
the so-called Zeno regimeii2iL& where the mixing time is 
much larger than the measurement (localization) time. 
In this regime, a good quantum readout measurement 
for an initial qubit state in the charge-state basis can 
be performed by repeatedly measuring the number of 
electrons, m{t) — Nji{t), that have tunneled through 
the right junction into the drain of the SET in the 
same detection time period t. One can then use the 
measurement results to construct the probability dis- 
tribution P{m,t). In Ref. Q, the probability distribu- 
tion P{m,t) = Trqb[pAr(m, i) + pN+iim^t)] is obtained 
from the Fourier analysis of the partially reduced den- 
sity matrix, Eq. (|30|l . The result, obtained in this way, 
is plotted in solid line in Fig. |21 This probability dis- 
tribution P{m, t) splits into two and their weights cor- 
respond closely to the initial qubit diagonal elements of 



p"(0) = 0.25 and p°°{0) = 0.75. 



In the quantum trajectory approach, P{m, t) can be 
explicitly simulated through constructing the histogram 
of the accumulated number of electrons Nuc = X) dN^c 
up to time t for many realizations of the detection records 
(generated together with their corresponding quantum 
trajectories), and then normalizing the distribution to 
one. The simulation of the normalized histogram shown 
in Fig.[21using 2000 quantum trajectories and their corre- 
sponding detection records is already in good agreement 
with the plot in solid line. However, the possible individ- 
ual realizations of measurement records and their corre- 
sponding quantum trajectories [e,g., in Fig.^ do provide 
insight into, and aid in the interpretation of the ensemble 
average properties. This is one of the appealing features 
of the quantum trajectory approach. 

For a charge qubit measured by a low-transparency 
point contact detector, this appealing feature of the 
quantum trajectories is illustrated in Ref. 8. Another ad- 
vantage of the quantum trajectory approach (or Bayesian 
formalism) is that it may describe a quantum feedback 
process. It has been shown^^ that one may utilize the 
measurement output for the feedback control and ma- 
nipulation of a qubit state. 



VII. CONCLUSION 

To summarize, we have derived the stochastic mas- 
ter equation for the SCB/SET system, which can be re- 
garded as a Monte Carlo method that allows us to simu- 
late the continuous quantum measurement process of the 
SCB qubit by the SET. We have shown that by taking 
a "partial" average over the fine grained measurement 
records of the tunneling events in the SET, this stochas- 
tic master equation reduces to the master equation pre- 
sented in Refs.[l|and|2. We have also presented numerical 
simulation for the dynamics of the qubit in a particular 
realization of the readout measurement experiment. We 
have shown that the probability distribution P{m, t) con- 
structed from 2000 quantum trajectories and their cor- 
responding detection records, is, as expected, in good 
agreement with that obtained from the Fourier analysis 
of the "partially" reduced density matrix. 
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